On the Dynamics of Unstable Quark-Gluon Plasma* 

Stanislaw Mrowczynski 

Institute of Physics, Jan Kochanowski University, 
ul. Swigtokrzyska 15, 25-406 Kielce, Poland 
and Soltan Institute for Nuclear Studies, 
ul. Hoza 69, 00-681 Warsaw, Poland 

(Received October 30, 209) 

Since the quark-gluon plasma, which is unstable due to anisotropic 
momentum distribution, evolves fast in time, plasma's characteristics have 
to be studied as initial value problems. The chromodynamic fluctuations 
and the momentum broadening of a fast parton traversing the plasma are 
discussed here. The two quantities are shown to exponentially grow in time. 
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1. Introduction 

The quark-gluon plasma (QGP), which is produced at the early stage of 
relativistic heavy- ion collisions, is most probably unstable due to anisotropic 
momentum distribution of quarks and gluons (partons), see the review fTj. 
The instability makes the system's state strongly time dependent, and con- 
sequently various plasma characteristics must be found as solutions of initial 
value problem. Two such characteristics are discussed in my lecture. I start 
with the fluctuations of chromodynamic fields and then, the correlation func- 
tions of the fields are used to compute the momentum broadening of a fast 
parton traversing the plasma. The two characteristics exponentially grow 
in time due to unstable modes. Throughout my lecture, which is based on 
two recent publications [21 [3], the plasma is assumed to be weakly coupled. 



2. Chromodynamic fluctuations 

In the quark-gluon plasma (QGP), which is on average locally color- 
less, chromodynamic fields, color charges and currents experience random 
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fluctuations. In the equilibrium plasmas there are characteristic stationary 
spectra of fluctuations which can be found by means of the fluctuation- 
dissipation relations. Fluctuations in nonequilibrium systems evolve in time, 
and their characteristics usually depend on an initial state of the system. 

Fluctuations can be theoretically studied by means of several methods 
reviewed in the classical monographs [H[5]. The method, which was chosen 
in [2], is clearly exposed in the handbook [6]. It is applicable to both equi- 
librium and nonequilibrium plasmas but the initial plasma state is assumed 
to be on average charge neutral, stationary and homogeneous. It will be 
shown that when the plasma state is stable, the initial fluctuations expo- 
nentially decay and in the long time limit one finds a stationary spectrum 
of the fluctuations. When the initial state is unstable, the memory of initial 
fluctuations is not lost, as the unstable modes, which are usually present in 
the initial fluctuation spectrum, exponentially grow. 

The chromodynamic fluctuations are studied here using the transport 
theory of weakly coupled quark-gluon plasma which is formulated in terms 
of particles and classical fields. The particles - quarks, antiquarks and gluons 
- should be understood as sufficiently hard quasiparticle excitations of QCD 
quantum fields while the classical fields are highly populated soft gluonic 
modes. The transport equation of quarks reads 



where Q(t, r,p) is the on-mass-shell quark distribution function which is 
N c x N c hermitean matrices belonging to the fundamental representation of 
the SU(iV c ) group; the covariant derivative in the four-vector notation reads 
D* 1 = — ig[A fl (x), ■ ■ ■ ] and E(i, r) and B(i, r) are the chromoelectric and 
chromomagnetic fields. The symbol {...,...} denotes the anticommutator. 
Since the fluctuations of interest are assumed to be of the time scale, which 
is much shorter than that of inter-parton collisions, the collision terms are 
absent in Eq. (TjQ). There are analogous transport equations for antiquark 
(<5(i,r,p)) and gluon (G(i,r,p)) distribution functions. 

The transport equations are supplemented by the Yang-Mills equations 
describing a self-consistent generation of the chromoelectric and chromo- 
magnetic fields by the color four-current = (p,j) 



where r a , T a with a = 1, ... , N~ — 1 are the SU(iV c ) group generators in the 
fundamental and adjoint representations. 

I consider small deviations from a stationary homogeneous state which 
is globally and locally colorless; there are no currents as well. The quark 
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distribution function of this state is Q^ m (p) = n(p)5 nm . Due to the absence 
of color charges and currents in the stationary and homogeneous state, the 
chromoelectric E(t, r) and chromomagnetic B(i, r) fields are expected to 
vanish while the potentials A°(t, r), A(i, r) are of pure gauge only. Since 
the plasma under considerations is assumed to be weakly coupled with the 
perturbative vacuum state, the potentials can be gauged away to vanish. 

The quark distribution function is written down as Q(t, r, p) = Q°(p) + 
5Q(t, r,p), and we assume that ^> \5Q\ and |V p (5°| 3> \V P 5Q\ with the 
analogous formulas for antiquarks and gluons. The transport ([1]) and Yang- 
Mills equations are linearized in deviations from the stationary homogeneous 
state. We assume that 5Q, E, B, ^4° and A are all of the same order. The 
linearized transport equation is 

d 

(di + v ' v ) 6Q ^' r ' p) ~ 9(E(t ' r) + v x B(t ' r)) Vpn(p) = ' 

After the linearization the Yang-Mills equations get the familiar form of 
Maxwell equations of multi-component electrodynamics. 

The linearized transport and Maxwell equations are solved with the 
initial conditions 5Q(t = 0, r,p) = SQo(r,p), E(i = 0,r) = Eo(r) and 
B(i = 0, r) = Bo(r), by means of the one-sided Fourier transformation 
defined as 

/*oo r 

/(u,,k) = J dt J d 3 re^- k - r )/(t,r). 

The chromoelectric field, which solves the equations, is found as 

[ - k 2 5 ij + WW + uj 2 e ij {uj, k)]Ei(u, k) = (2) 

_ f d 3 p ^(v x B a0 (k)) J V^(p) 
% 2 J (2vr) 3 u - v • k 

- 9" J ^ k ,v 6N $ ( k > P) + ^o(k) - *(k x B a0 (k)Y , 

where f(p) = n(p) + n(p) + 2N c n g (p) and <5iV a (r,p) = Tr[r a (<5Q (r, p) - 
5Qo(r, p)) +T a 5Go(r, p)]; e J:, (u;,k) is the chromodielectric tensor of, in gen- 
eral, anisotropic plasma in the collisionless limit; e l:, (u;,k) does not carry 
any color indices, as it corresponds to a colorless state of the plasma. 

When the plasma stationary state is isotropic, the dielectric tensor can 
be expressed through its longitudinal (ex(cj, k)) and transverse (er^k)) 
components and the matrix £ lJ (u;,k) = — k 2 ^ + WW + a; 2 e*- ? (a;, k) from 
the left-hand-side of Eq. ([2]) can be inverted as 

<s->)<V,k) = -^-^ + , , » V M > - . 

o^et^o^k) k z uj z st{uj, k) — k^ V k z / 
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Then, Eq. ([2]) provides an explicit expression of the chromoelectric field. 
Using the Maxwell equations, the chromomagnetic field, color current and 
color density can be all expressed through the chromoelectric field. 

The correlation functions (E l a (ti, r{)El (t 2 , r 2 )}, {B l a (ti, vi)B 3 b {t 2 , r 2 )), 
where (• • •) denotes averaging over statistical ensemble, are determined by 
the initial correlations such as (SNg (r±, pi)5N^ (r 2 , p 2 )), (Ea (ri)E J bQ (r 2 )), 

(5Nq (n, pi)£ , ^ (r 2 )} which can be all expressed, using the Maxwell equa- 
tions, through the correlation function of the distribution functions. The 
latter one is identified with the respective correlation function of the classical 
system of free quarks, antiquaries and gluons which on average is stationary 
and homogeneous. For quarks the free correlation function is 

(5Q mn (t 1 ,r 1 ,p 1 )5Q pr (t 2 ,r 2 , P2 )) hec = 5 mr ^(2vr) 3 ^( Pl - p 2 ) 

x ^ 3 )(r 2 - n - vi (i 2 - h)) n(pi) . 

In the case of equilibrium plasma, where all collective modes are damped, 
I consider the times which are much longer than the decay time of collec- 
tive excitations. Then, the correlation function of the chromoelectric fields 
equals 

(^( il ,r 1 )^(t 2 ,r 2 )) 0O = J ^^03e-(^^)~ k -(---))(^^') fc , 
where the fluctuation spectrum is 



(Ei4) k = 25 ab Tco 3 



+ [0 J 



k 2 |c; 2 e L (cj,k)| 2 V k 2 J |u; 2 e T (u;,k) -k 2 | 2 



As seen, the fluctuation spectrum has strong peaks corresponding to the 
collective modes determined by the equations £c(o>, k) = and u> 2 eT(uj, k) — 
k 2 = 0. 

As an example of a nonequilibrium situation, I discuss fluctuations of 
longitudinal chromoelectric fields in the two-stream system which is unstable 
with respect to longitudinal modes. Nonequlibrium calculations are much 
more difficult than the equilibrium ones. The first problem is to invert the 
matrix T, ij (lo, k). In the case of longitudinal electric field, which is discussed 
here, the matrix is replaced by the scalar function. 

The distribution function of the two-stream system is chosen to be 

/( P ) = (2^) 3 nr5( 3 )( P -q)+5( 3 )(p + q) 



where n is the effective parton density in a single stream. There are four 
roots ±w±(k) of the dispersion equation ei(tj,k) = 0. The solution u;+(k) 
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represents the stable modes and w_(k) corresponds to the well-known two- 
stream electrostatic instability for ku ^ and k 2 (k-u) 2 < 2/x 2 (k 2 — (k-u) 2 ) 
where u = q/|q| is the stream velocity and /j, 2 = g 2 n/2\q\. Then, cj_(k) = 
*7k with < 7k € R. 

The correlation function of longitudinal chromoelectric fields generated 
by the unstable modes is found as 

rO^^ra))^^ = (3) 

(? 6 ab n f d3k e ik ( ri - r2 )FF( 7 2 + (k-u) 2 ) 2 



2 J (2vr)3 k 4 {uo\ - u?_ ) 2 7 2 

( 7 2 + (k • u) 2 ) cosh ( 7k (t! + t 2 )) + ( 7 2 - (k • u) 2 ) cosh ( 7k (t! - t 2 )) 

As seen, the correlation function of the unstable system is invariant with 
respect to space translations - it depends on the difference (ri — Y2) only. 
The plasma state, which is initially on average homogeneous, remains like 
this in course of the system's temporal evolution. The time dependence of 
the correlation function is very different from the space dependence. The 
electric fields exponentially grow and so does the correlation function both 
in (t\ + t 2 ) and (t\ — t 2 ). The fluctuation spectrum also evolves in time, 
as the growth rate of unstable modes is wave-vector dependent and after a 
sufficiently long time the fluctuation spectrum is dominated by the fastest 
growing modes. 



3. Momentum broadening of a fast parton 

When a highly energetic parton travels through dense QCD matter, 
it receives random kicks from elastic interactions with constituents of the 
plasma. The average transverse momentum transfer per unit path length 
is related to the radiative energy loss of the parton [7\. The parameter 
describing the average amount of transverse momentum broadening per unit 
length is called q and is defined as 

q = d(Ap 2 T )/dz, (4) 

when the fast parton flies along the direction z. The values of q extracted 
from experimental data on relativistic heavy-ion collisions vary in a rather 
broad range 0.5 — 15 GeV 2 /fm depending on the model of hard particle 
propagation in strongly interacting matter produced in nuclear collisions 

mm- 

The calculations of q for the case of perturbative quark-gluon plasma in 
equilibrium are well understood (see [10\ PTT] for recent work) . For such a 
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plasma the value of q is predicted to lie at the lower end of the range of values 
deduced from experiments [9]. However, the plasma momentum distribution 
is initially anisotropic and recently q has been computed |121 [13] for such 
a plasma. However, the fact that the anisotropic plasma as an unstable 
system evolves fast in time has not been taken into account. It seems 
rather unjustified, as the numerical simulations |14[|15j clearly indicate that 
q receives a sizable contribution from the unstable growing modes. 

An analytic approach to compute q in unstable plasma has been devel- 
oped in [3]. It formulates the transverse momentum fluctuations in terms 
of classical Langevin problem, q is computed by treating the parton as an 
energetic classical particle with SU(3) color charge moving in the presence 
of the fluctuating color fields. Then, q is expressed through the correla- 
tion function of chromodynamic fields computed in [2] . For the equilibrium 
plasma the Langevin approach recovers the known result obtained within 
the standard thermal field theory |16j . 

Let me consider a classical parton which moves across a quark-gluon 
plasma. Its motion is described by the Wong equations [T7] 



^ = (5) 



(It 



gQ a {T)Fr(x{T))u u {rl (6) 



= - 9 r bc P,(r)A^x(r))Q c (r), (7) 



where r, x^(t), u^(t) and p^{r) are, respectively, the parton's proper time, 
its trajectory, four-velocity and four-momentum; Fa and A% denote the 
chromodynamic field strength tensor and four-potential, respectively, and 
Q a is the classical color charge of the parton. 

We look for a solution of the Wong equations in a specific gauge assuming 
that the potential vanishes along the parton's trajectory i.e. our gauge 
condition is 

p,(r)A^x(r)) = 0. 

Then, Eq. <^ simply tells that Q a is constant as a function of t. 

One solves Eqs. (0EI) assuming that the parton's momentum p is so high 
that its changes Ap caused by the interactions with the medium are small 
compared with p. The changes of the velocity vector v are then negligible, 
and we can consider the parton to move along a straight-line path with 
constant velocity. 

Assuming that the parton moves with the speed of light in the positive 
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z-direction that is x^(t) = (i, 0,0,i), one finds 

pP(jt) = p«(0) + gQ a I dt' [Ff(t>) - Ff(t)] , (8) 
Jo 

where Fa h (t) should be understood as a short hand notation of Fa U (x(t)). 

In the spirit of Langevin approach, we consider the ensemble average 
{p^ '{t)p h '{t)) indicated by the angular brackets. The ensemble average in- 
volves averaging over color charges which is performed by means of the 
relation 

' dQQ a Q b = C 2 5 ab 



where C 2 = 1/2 for particles (quarks) in fundamental representation of the 
SU(iV c ) group and C 2 = iV c for particles (gluons) in adjoint representation. 
Then, we find 



(PW(*)> = <P"(0)p"(0)> (9) 
+9 2 J^J*dtiJ*dt 2 {(F^(t 1 )-Ff(t 1 ))(F^ (t 2 )-Ff (i 2 ))} , 

where Cr (R = F, A) is the eigenvalue of the quadratic Casimir operator, 
C F = (7V C 2 - l)/2N c and C A = N c . 

v. Ml .,,„•! / \r,. ( 



Introducing Ap T (i) = Pt(*)~Pt(0) and (Apl(t)) = (Ap T (t)-Ap T (t)), 



we have 



- {EZ(h)BV(t 2 )) + (Ey(h)B-(t 2 )) - (BV^E^t,)) (10) 
+ (B*(h)Ey(t 2 )) + (B* a (h)B* a (t2)) + (Bl^Bl^)) 

where, say, E£(t) should be understood as E%(t,r(t)) with r(t) = (0, 0, t) 
that is only the fields at the parton's trajectory enter Eq. (|10p . 

Eq. PU|) has been derived in the specific gauge and the right-hand side of 
the equation is, in general, gauge dependent. However, the field correlation 
functions derived in [2] are gauge independent within the Hard-Loop Ap- 
proximation, as discussed in detail in Sec. VIIIA of [2]. Therefore, Eq. (fTUj) 
is gauge independent within the used approximations. 

Let us now assume that the quark-gluon plasma is translationally in- 
variant in space and time and, hence, the field correlators depend only on 
the difference of the field's arguments. This assumption, which is relevant 
for equilibrium plasmas, will be not adopted for the two-stream system as 
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the growth of unstable modes break the translational invariance in time. 
Making use of the translational invariance, one introduces the fluctuation 
spectrum 



(KEi) k = JdtJ d 3 re 4 ^~ kr )(^(t,r)^(0,0)} 



11) 



which allows us to write down Eq. (fTOj) in the form 
<Ap^)> = g 2 J^-r tdh fdt 1 ' 1>h ' 



N?-l 



o 



o 



(2k)* 



{E x a E%) k + (Ej>El>) h - (E*B») k + {Ej>I%) k 
■(B»E%) k + {BlEl) k + (B* a B* a ) k + {B y a B y a )y 



(12) 



Since the double integral over t\ and t 2 tends to a delta function of (u — k z )/2 
in the long-time limit 



4 sin 



(w — k z )t 



dh / dt 2 e i ^- k ^ tl ' t ^ 
ii Jo (w - k z y 



'— ^ L?( u /••• 



we find the transport coefficient q = d{Apj,(t))/dt as 



q = 2itg 



Cr 



d A k 



8(u) - k z ) 



{E*E* a ) k + {ElEl) k 



N*-lJ (2tt)4 

(EZBl) k + (i^)* - {BlE x a ) k + (S^) fc + + (fl^)* 



(13) 



Using the equilibrium field correlators derived in j2J 3 Eq. (|13|) gives 



9 = 2g 2 C R T 



d 3 k 



(2vr) 3 ^k 2 



%e L (k z ,k) k 2 z k 2 %e T (k z ,k) 



+ 



£ L (A; 2 ,k)| 2 |fcf £T (fc Z) k)-k 2 | 2 



(14) 



The classical formula (|14p holds for the inverse wave vectors of the fields 
which are much longer than the de Broglie wavelength of plasma particles. 
It requires |k| <C T where T is the plasma temperature. For larger wave 
vectors a quantum approach discussed in [3] is needed. 

Let us now consider the momentum broadening of a fast parton in un- 
stable anisotropic plasmas. For the sake of analytical tractability, the two- 
stream plasma is discussed and only longitudinal electric fields are taken 
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into account. Substituting the correlation function ([3]) into Eq. (I10p . one 
finds 



9 



C R n 



d 3 k 



dt 7 



krp 



(^(ti) El{t 2 )) + {El{t x ) El{t 2 )) 
ill + (k • u) 2 ) 2 



(2vr) 3 k 4 (o;2 



2 \2 



(15) 



( 7 £ + (k-u) 2 )(|e^ +7k )' 
7k 



1| 2 _|_ | e (* fc z-7k)* _ ]|2\ 



+4( 7 ^ - (k ■ u)') 



k 2 - ~2 



2n "'z 



^ 2 +7k. 

When only the fastest growing mode is included, Eq. (|15p changes to 



rft ~ 2 



(2^)3 



2\3 



4( 7 2 + (k • u) 2 ) 



2 \2 



7k(fc 2 + 7k) ' 



(16) 



It should be noted that Eq. (fT6j) . in contrast to Eq. (|T5|) . suffers from diver- 
gence when 7 k — > 0. The momentum broadening (|16p grows exponentially 
in time, as the exponentially growing fields exert an exponentially growing 
influence on the propagating parton. This effect is missing in the previous 
results [121 13] obtained for the anisotropic plasma treated as a stationary 
medium. 



4. Discussion and Outlook 

If the momentum distribution of partons is anisotropic, the quark-gluon 
is unstable with respect the chromomagnetic modes. The instability growth 
rate 7 is of order gT [1J where T is here the characteristic parton momentum 
corresponding to temperature of equilibrium plasma. The time scale of 
unstable mode growth 7 -1 is the shortest dynamical time scale in weakly 
coupled plasmas. Therefore, the unstable plasma cannot be treated as a 
static medium whenever a characteristics, which involves color degrees of 
freedom, is studied. 

In this lecture I discussed the fluctuations of chromodynamic fields in 
the plasma, and then, the field correlation functions were used to compute 
the momentum broadening of a fast parton flying across the plasma. Both 
quantities were found as solutions of initial value problem. In the case of 
unstable plasma, the field fluctuations as well as the momentum broadening 
were explicitly shown to exponentially grow in time. 
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An analytic treatment of unstable plasma is rather difficult. Therefore, 
instead of the plasma with momentum distribution relevant for relativistic 
heavy-ion collisions, a toy model representing the two stream system was 
discussed. The model grasps some important features of unstable systems, 
but it cannot be used for any quantitative predictions on the quark-gluon 
plasma produced in high-energy nucleus-nucleus collisions. An analysis the 
chromodynamic fluctuations and momentum broadening in such plasma is 
in progress but many other plasma characteristics, in particular transport 
coefficients, need to be studied to understand QGP from early stage of 
relativistic heavy-ion collisions. 
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